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Abstract 

This paper is concerned with a fluidodynamic model for traffic flow. More precisely, we 
consider a single conservation law, deduced from conservation of the number of cars, 
defined on a road network that is a collection of roads with junctions. The evolution 
problem is underdetermined at junctions, hence we choose to have some fixed rules for 
the distribution of traffic plus an optimization criteria for the flux. We prove existence, 
uniqueness and stability of solutions to the Cauchy problem. 

Our method is based on wave front tracking approach, see and works also for 
boundary data and time dependent coefficients of traffic distribution at junctions, so 
including traffic lights. 
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1 Introduction 



This paper deals with a fluidodynamic model of heavy traffic on a road network. More 
precisely, we consider the conservation law formulation proposed by Lighthill and Whitham 
[]l2| and Richards [|l3) . This nonlinear framework is based simply on the conservation of cars 
and is described by the equation: 

Pt + f( P ) x = 0, (f.l) 

where p — p(t,x) £ [0, p ma x}, (t,x) £ M 2 , is the density of cars, v{t,x) is the velocity and 
f{p) = v p is the flux. This model is appropriate to reveal shocks formation as it is natural 
for conservation laws, whose solutions may develop discontinuities in finite time even for 
smooth initial data, [||. In most cases one assumes that v is a function of p only and that 
the corresponding flux is a concave function. We make the same assumption, moreover we 
let / have a unique maximum a £]0, p m ax[ and for notational simplicity assume p ma x = I- 

Here we deal with a network of roads, as in |lf} |. This means that we have a finite number 
of roads modeled by intervals [ai,&i] (with one of the two endpoints possibly infinite) that 
meet at some junctions. For endpoints that do not touch a junction (and are not infinite), 
we assume to have a given boundary data and solve the corresponding boundary problem, 
as in H ||, ||, |H . The key role is played by junctions at which the system is underdetermined 
even after prescribing the conservation of cars, that can be written as the Rankine-Hugoniot 
relation: 

n n+m 
i— 1 j—n+1 

where pi, i = 1, . . . , n, are the car densities on incoming roads, while pj, j = n + 1, . . . ,n + m, 
are the car densities on outgoing roads. In fl"lf| , the Riemann problem, that is the problem 
with constant initial data on each road, is solved maximizing a concave function of the fluxes 
and it is proved existence of weak solutions for Cauchy problems with suitable initial data of 
bounded variation. In this paper we assume that: 

(A) there are some prescribed preferences of drivers, that is the traffic from incoming roads 
is distributed on outgoing roads according to fixed coefficients; 

(B) respecting (A), drivers choose so as to maximize fluxes. 
To deal with rule (A), we fix a matrix 

X^jiJj— ...n+m, i— l,...,n ^ ^' , 

such that 

n+m 

exji^aji', < ctjn < I, 2J a ji = !> (!- 3 ) 

j=n+l 



I 
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for each i' ^ i = 1, ... , n and j = n+ 1, n+m, where Ojj is the percentage of drivers arriving 
from the i— th incoming road that take the j — th outgoing road. Notice that with only the 
rule (A) Riemann problems are still underdetermincd. This choice represents a situation 
in which drivers have a final destination, hence distribute on outgoing roads according to a 
fixed law, but maximize the flux whenever possible. We are able to solve uniquely Riemann 
problems and, in case of simple junctions with two incoming and outgoing roads, to generate 
a Lipschitz semigroup, defined on L , whose trajectories are weak solutions and respect rules 
(A) and (B) in case of bounded variation (the same conditions are not meaningful if the 
solution is only L 1 ). Our main technique is the use of a front tracking algorithm and suitable 
approximations and functionals to control the total variation. We refer the reader to || for 
the general theory of conservation laws and for a discussion of wave front tracking algorithms. 

The main difficulty in solving systems of conservation laws is the control of the total 
variation, see [0]. It is easy to see that for a single conservation law the total variation is 
decreasing, however in our case it may increase due to interaction of waves with junctions. 
The problem is quite delicate, as shown in Appendix B, where an example is given of a single 
wave of arbitrarily small strength (variation) that, interacting with a junction, generates 
waves whose strengths are bounded away from zero. Hence we can not expect any bound on 
the total variation of the solution in term of the total variation of the initial data, as it is 
the case for systems. This arbitrarily large magnification of total variation is possible only if 
waves crossing the value a interact with junctions at which the boundary data of the roads 
are bad, that is in [0, a] for incoming roads and in [a, 1] for outgoing roads. We thus have 
first to deal with special data of bounded variation that have a finite number of crossing of 
the value a. The sum of the number of these crossing plus the number of bad boundary data 
is proved to be decreasing along front tracking approximate solutions. 

However, this is not enough since the variation can still increase due to interactions with 
junctions (and there is no bound on the number of interactions). The conserved quantity 
is the total variation of the flux. We prove this fact for junctions with only two incoming 
roads and two outgoing ones, and show in Appendix A that the conclusion does not hold 
for for junctions with three (or more) incoming and outgoing roads. Unfortunately the total 
variation of the flux is not equivalent to the total variation of p, since f'(a) = 0. We thus 
have to approximate the flux with one having never vanishing derivative and a corner at a, 
and then pass to the limit. 

Our techniques are quite flexible, so we can deal with time dependent coefficients for 
the rule (A). In particular we can model traffic lights and also in this case the control of 
total variation is extremely delicate. An arbitrarily small change in the coefficients can 
produce waves whose strength is bounded away from zero. Still it is possible to consider 
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periodic coefficients, a case of particular interest for applications. We can also deal with 
roads with different fluxes: this can be treated in the same way with the necessary notational 
modifications. 

There is an interesting ongoing discussion on hydrodynamic modelization for heavy traffic 
flow. In particular some models using systems of two conservation laws have been proposed, 
see (| !§, 0. We do no treat this aspect. 

The paper is organized as follows, in Section |2] we give the definition of weak entrtopic 
solution and following (A) and (B) we introduce an admissibility condition. In Section [| we 
prove the existence and uniqueness of admissible solutions for the Riemann Problem in a 
junction, then using this we describe the construction of the approximants for the Cauchy 
Problem (see Section [|) . In section |^ we prove the monotonicity of the number of big waves 
for piecewise constant solutions. Assuming that /' is bounded away from and that there 
are at most two incoming and outcoming roads in each junction we prove the monotonicity 
of the total variation of the flux (see Section |^) and existence, uniqueness and stability of 
admissible solutions for the Cauchy Problem with suitable BV initial data. Using these 
results we show the existence of a unique Lipschitz semigroup of solutions defined on L 1 (see 
Section |J) also in the case in which / is smooth. In Section ^ we describe what happens 
when there are traffic lights and time dependent coefficients. In Appendix ^ we show with an 
example that the total variation of the flux can increase when there are three incoming and 
three outcoming roads in a junction. Finally, in Appendix [b| we show that the interaction of 
a small wave with a junction can produce a uniformly big wave. 



2 Basic Definitions 

We consider a network of roads, that is a modeled by a finite collection of intervals U = 
[ai, bi] C M, i = l,..., JV, possibly with either a t — — oo or b t — +oo, on which we consider 



the equation (1.1). Hence the datum is given by a finite collection of functions pi defined on 
[0, +oo[xIj. 

On each road Ii we want pi to be a weak entropic solution, that is for <p : Zj — > 1R smooth 
with compact support on ]0, +oo[x]aj, bi[ 



(pi^ + f( P i)^)dxdt = 0, (2.4) 



+ °° I** ( dip ,. 
/ \P 

and for every k S M and (p : ij —> 1R smooth, positive with compact support on ]0, +oo[x]<2j, 6, 



o 



(Vi-fel-^+sgn { Pi -k)(f( Pi )-f(k))-£)dxdt>0, (2.5) 
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It is well known that for any initial data in L°°, denned on the whole M, there exists a 
unique weak entropic solution depending in a Lipschitz continuous way from the initial data 
in L\ oc . 

We assume that the roads are connected by some junctions. Each junction J is given by 
a finite number of incoming roads and a finite number of outgoing roads, thus we identify 
J with ((ixj • • • ; i n )i (ii) • • • ijm)) where the first n-tuple indicates the set of incoming roads 
and the second m-tuple indicates the set of outgoing roads. We assume that each road can 
be incoming road at most for one junction and outgoing at most for one junction. 

Hence the complete model is given by a couple (I, J), where X = {Jj : i = 1, . . . , N} is 
the collection of roads and J the collection of junctions. 

Fix a junction J with incoming roads, say I\, . . . , J„, and outgoing roads, say I n +i , ■ ■ ■ , I n +m ■ 
A weak solution at the junction J is a collection of functions pi : [0, +oo[x/ i — > ]R, i — 
1, . . . , n + m, such that 

for each ipi, ip n +m smooth having compact support in ]0, +oo[x IR, that are also smooth 
across the junction, i.e. 

d^p dtp 



Remark 2.1 Let p = (p±, . . . , p n + m ) be a solution of (LI) such that each x — * pi(t, x) has 
bounded variation. We can deduce that it satisfies the Rankine-Hugoniot Condition in the 
junction J, namely 

n n+m 

£/(*(*,&<))= E /(Pi(*.Oj))> (2-7) 

i— 1 j=n+l 

for each t > 0. 



Remark 2.2 The assumption 7^ a^j/ in (L3)is needed in order to have uniqueness of so- 



lutions to Riemann problems, see Section [|. However, this condition can be relaxed requiring, 
for example, that if ctji = atji' then the fluxes from road i and i' coincide. 

The rules (A) and (B) can be given explicitly only for solutions with bounded variation 
at each time as in next definition. 



Definition 2.1 Let p = (pi, . . . , p n +m) be such that pt(t, ■) is of bounded variation. Then p 



is an admissible weak solution of (1.1) related to the matrix A satisfying (l.c) at the junction 
J if and only if following properties hold: 
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(i) p is a weak solution at the junction; 

n 

fti) f(Pj(-> a 3+)) = J2 otjif(Pi{',bi-)), for each j = n+ 1, ...,n + m; 

i=l 

n n+m 

(Hi) f(fli(-,bi—)) + f(pj(-,aj+)) is maximum subject to (ii). 

i—l j—n+1 

For every road Ii = [oj, bi], if a, > — oo and Ii is not the outgoing of any junction, or b t < +00 
and Ii is not the incoming road of any function, then a boundary data ipi : [0, +00 [— > JR is 
given. In this case we ask pi to satisfy pi(t,di) — ipi(t) (or pi(t,bi) = ipi(t)) in the sense of 
I- 

Our aim is to solve the Cauchy problem on [0, +00 [ for a given initial and boundary data 
as in next definition. 

Definition 2.2 Given pi : Ii — > M and possibly ipi : [0,+oo[^ M, functions of L°° , a 
collection of functions p — (p l5 . . . , Pn) with pi : [0, +oo[xi i — > 1R continuous as functions 



from [0,+oo[ into Lj oc , is an admissible solution if Pi is a weak entropic solution to (1.1) on 
Ii, pi(0,x) = pi{x) a.e., pi(t,bi) — ipi{t) ^ n the sense of Mj ; finally such that at each junction 
p is a weak solution and is an admissible weak solution in case of bounded variation. 

The treatment of boundary data in the sense of [|| can be done in the same way as in |2[ || , 
thus we treat the case without boundary data. All the stated results hold also for the case 
with boundary data with the obvious modifications. 
On the flux / we make the following assumption 

(T) f : [0, 1] — * M is smooth, strictly concave (i.e. /" < — c < for some c > 0), 
/(0) = /(l) = 0, |/'(a;)| < C < +00 and there exists a e]0, 1[ such that f (a) = (that 
is a is a strict maximum). 



3 The Riemann Problem 

In this section we study Riemann problems. For a scalar conservation law a Riemann problem 
is a Cauchy problem for an initial data of Heaviside type, that is piecewise constant with 
only one discontinuity. The solutions of these problems are the building blocks to construct 
solutions to the Cauchy problem via wave front tracking. These solutions are formed by 
continuous waves called rarefactions and traveling discontinuities called shocks. The speed 
of waves arc related to the values of f'(p). 

Analogously, we call Riemann problem for the road network the Cauchy problem corre- 
sponding to an initial data that is piecewise constant on each road. The solutions on each 
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road Ii can be constructed in the same way as for the scalar conservation law, hence it suffices 
to describe the solution at junctions. Because of finite propagation speed, it is enough to 
study the Ricmann Problem for a single junction. 

As explained in the Introduction, we first have to treat the case of fluxes with nonvanishing 
derivative, hence we assume that 

(JF1) / : [0, 1] — ► M is continuous, strictly concave, /(0) = /(f) = and there exists a e]0, 1[ 
such that / is smooth on [0, a[ and on }a, 1] and 

< c < \f'(x)\ < C < +oo, (3.8) 

for each x G [0, er[U]cr, 1]. 

Consider a junction J in which there are n roads with incoming traffic and m roads with 
outgoing traffic. For simplicity we indicate by 

(t,x) e M+ x Ii t->pi(t,x) e [0, 1], i = l,...,n (3.9) 

the densities of the cars on the road with incoming traffic and 

(t, x) € JR+ x Ij i— ► Pj(t, x) € [0, 1], j = n+ 1, ...,n + m (3.10) 

those on the roads with outgoing traffic. 




figure f 

We need some more notation: 
Definition 3.1 Let t : [0, f] — > [0, 1], t(<t) = a, be the map satisfying the following 

r{p) * p, f(r(p)) - f(p), 

for each p ^ a \ 

Clearly r is well defined and satisfies 

< p < a t=> a < t( P ) < 1, a < p < f < r(p) < a. 
The main result of this section is the following Theorem. 
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Theorem 3.1 Let f : [0, 1] — ► 1R satisfy (J-l) and p\$, p n +m,o G [0, 1] be constants. 
There exists ctn unique admissible weak solution, in the sense of D efinition \Z. A p — (pi? •••) Pn-\-r, 



of (1.1) at the junction J such that 

Pl(0, ') = Pl,0, , Pn+m(0, •) = Pn+m.O- 

Moreover, there exist a unique {n + m) — tuple (p%, ...,p n+m ) G [0, l] n + m such that 



Pi G 



{Pi,o}U]r(pi i0 )l] if0< p lfl < a, 
[o-, 1] if o- < p ifi < 1, 



i = 1, n 



and 



Pj G 



[0,a] ifO< Pj ,o<(T, 
{Pj,o} U [0, r(p j)0 )[ «/cr < Pj,o < 1, 



j = n + 1, n + m. 



Fixed i G {1, ...,«}, if Pi,o < /3j t/iere results 

( p. Q lfx < /fe)-/(p,o) 



a;) 



Pi-Pi,0 



Pi otherwise, 



and if p. i < 

!Pi,0 if x < f'(pifl)t, 

(f)~\f) iff(pi,o)t<x<f'(pi)t, 

Pi ifx>f'(pi)t. 

Otherwise, fixed j G {n + 1, n + m}, if pj$ < /5j </iere results 

Pj if x < f'{p.j)t, 

Pi&x) = \ {f'T\l) iffiPo)t < * < f(pj,o)t, 

Pj.o ifx>f'(pj, )t, 



and if pj < pj.o 



Proof. Define the map 



Pj(t, x) 



[Hi ifx< f ^_-^\ 
y Pjfi otherwise. 



and the sets 



E:( 11 ,..., ln )eM n ^J2^ 

i=l 

[0, f{p i>0 )] if < p lfi < o-, 
[0,/(a)] ifo-<p it0 <l, 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
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n = 



. j [0,f(*)} ifO< ft , <a, . 
\L ■ = < 7 = n + 1, n + m, 

\ [QJ( Pj , )] if*< Pj , <l, ■ 

{(7i, -)7n) S Oi x .... x Sl n \A ■ (71, ...,7„) T € 0„+i x .... x fi n+m j. 



Since -E is linear, the set f2 is closed, convex and not empty. By (1.3) there exists a unique 
vector (71, ...,7„) £ such that 

£(71, ...,7„) = max E(ii, ...,7„). 
(7iv,7«)en 

Fix i G {1, ...,n}, let pi £ [0, 1] be such that 

f(Pi) =% Pi e j 
By (^"1), /5i exists and is unique. Let 

n 

lj=22<xji% j = n + l,...,n + m 
and /5j G [0, 1] be such that 



{pi,o}U]T(p ifi )l} if < p i: o < cr, 
[(7,1] if cr < /jj.0 < 1. 



/(ft)=7j ; e 



[0, <t] if < Pj - o < <r, 

{p i: o} U [0,r(pj,o)[ if cr < Pjfi < 1. 

Since (71, 7„) € f2, pj exists and is unique. Solving the Riemann Problem (see || Chapter 
6]) on each road, the thesis is proved. □ 



4 The Wave Front Tracking Algorithm 

Once the solution to a Riemann problem is provided, we are able to construct piecewise 
constant approximations via wave-front tracking. The construction is very similar to that 
for scalar conservation law, see (H, hence we only briefly describe it. 

Let po be a piecewise constant map defined on the road network. We want to construct 



a solution of (1.1) with initial condition p(0, •) = po- We begin by solving the Riemann 



Problems on each road in correspondence of the jumps of po and the Riemann Problems at 



junctions determined by the values of po (see Theorem 3T). We split each rarefaction wave 
into a rarefaction fan formed by rarefaction shocks, that are discontinuities traveling with 
the Rankine-Hugoniot speed. We always split rarefaction waves inserting the value a (if it is 
in the range of the rarefaction), in order to control the number of big waves defined in next 
Section. 
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When a wave interacts with another one we simply solve the new Riemann Problem, if 
otherwise it reaches a junction then we solve the Riemann Problem at the junction. Since 
the wave speed is bounded there are finitely many waves on the network at each time t > 0. 
We call the obtained function a wave front tracking approximate solution. Given a general 
initial data, we approximate it by a sequence of piecewise constant functions and construct 
the corresponding approximate solutions. If they converge in Lj oc , then the limit is a weak 
entropic solution on each road, see M for a proof. 



5 Estimates on the Number of Big Waves 

In this Section we consider big waves, that are the waves crossing the value a. For these 
waves the variation of f(p) is not comparable to the variation in p, more precisely the former 
can vanish while the second is different from zero. Since only the variation of f{p) happens 
to be conserved we need to control the number of big waves. 

Let p a piecewise constant map defined on the network and J a junction with n roads 
with incoming traffic and m roads with outgoing traffic as in Section |^. Define the set 

$j(p) = {i e {l,...,n}\pi(bi-) G [0,cr]} U {j G {n+l,...,n + m}\pj(aj+) G [<r,l]}. 

For each road Ii we denote by {x a G h : a £ Ai} and by {(p(x a —),p(x a +)) : a £ Ai} the 
set of discontinuity points and the set of discontinuities, respectively, of the map pi on the 
road Ii. We define 

G l {p) = {(p(xa-),p(x a +))\ a G Ai, sgn (p k (x a -) - a) ■ sgn (p k (x a +) - a) < 0}, 
with the agreement sgn (0) = 0, and the functional 

N 

A^)=£#($ 7 (p))+]T#(G>)), 

JdJ i=l 

where # indicates the cardinality of a set. 

The main result of this section is the following. 



Lemma 5.1 Let f : [0, 1] — ► M satisfy (J-l) and p be a piecewise constant wave front tracking 
approximate solution of (1A) on the net. Then the map 

t >0i — > N(p(t,-)) 



does not increase. 
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Proof. We begin considering a single junction J with n roads with incoming traffic and 
m roads with outgoing traffic and an equilibrium configuration (pi,o, Pn+m,o) G [0, l]™+ m , 
namely the solution of the Riemann Problem in the junction with that initial data is constant. 
Suppose that a wave on one road arrives to the junction at time t and there is no other wave 
on the roads, then we claim that N(p(i—, •)) = N(p(i+, •)). Assume that the wave is on an 
incoming road, for example the first one and let (pi,pi,o) be the values on the left and right 
side of the wave respectively. Since the wave is approaching the junction, its speed is positive 
and so < p\ < a, moreover since it is the unique wave 

G*(p(F- •)) =0, i = 2,...,n + m. 

Let (pi, Pn+m) be the solution to the Riemann Problem with initial data (pi, P2,o---j Pn+m,o) 
(see Theorem 3.1), there results 

pi G { Pl } U [a, 1] 

Pi G {Pi,o} U [a, 1], i = 2, ...,n 

Pj £ {Pj,o} U [0,ct], j = n + 1, ...,n + m. 

In the following we study the change of the functional N due to the presence of new waves. 
If a new rarefaction is produced then it can not cross the value <r, otherwise there would 
be rarefaction shocks with positive and negative velocity at the same time. Hence each 
functional G l can not be bigger than one after the interaction. By abuse of notation, we 
indicate the whole rarefaction fan as a single wave for notational simplicity. 
So, fixed i G {2, n}, if % $j(p(i—, •)), then 

«^j(p(t+,-)) 1 (Pi,o,Pi) ?G*(p(i+,-)), 

and if i G &j(p(t— 7 •)) we have 

k = Pifi i G •)). (Pi,o,Pi) $ G l (p(t+, •)), 

o < Pi < 1 i & <t>.j{p(i+, ■)), (Pi,o,Pi) G G l (p(t+, •)). 
On the other hand, fixed j G {n+ 1, n + m}, if j £ , •)), then 

j ? $j(p(t+,-)), (pj,ao) 2 G?(p(t+,-)), 

and if j G $j( ( o(f— , •)) we have 

Pi = Pi,o =^J£ (p(<+i 0) , (Pj , ft.o) ^ •)) , 

< < a j ^ $./( j0 (t+, •)), (Pi,Pi,o) G G" (/»(*+, •))■ 
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Hence the contribution to N due to roads Jj, i = 2, . . . , n + m, does not increase. Let us now 
treat the waves on the first road. 

Notice that if pi_ = a then p\ ^ a hence 

1 G Mp(*-, ■)). (Pi,Pi,o) G G 1 (p(i-, •)) 

and iV can not increase. The same conclusion holds if < pi,o < u and p\ = a . 
Now, if < pifi < a and p\ ^ <r, then there results 

1 G ■)), (pi,Pi,o) £ G 1 (p(i-, •)) 

and 

Pi =Pl => 1 G $j(p(t+,-)), (Pl.Pl) ^ •)), 

Pi ^pi =*<7 < pi < 1 => 1 £$j(p(i+,-)), (Pl,/3l) G G 1 (p(f+,-))- 

If it < pi^o < 1 and pi = g we have 

1 £ Mp(*-, ■)). (Pi,Pi,o) G •)) 

and 

Pi =Pi =ct^ 1 G *j(p(t+,-)), (Pi)Pi) ^ •)), 
<r < Pi < 1 => 1 £ $j(p(i+,-)), (pi,Pi) G G 1 (p(i+, •))• 
Finally, if cr < pi ; o < 1 then p\ ^ a, pi ^ a, and wc have 

1 ^*j(p(t- ■)), (pi,Pi,o) G G 1 (p(i- 7 •)), 

hence 

Pi =Pi =► 1 G (pi.Pi) 2G 1 (p(t+,-)), 

<r < Pi < 1 1 £ *j(p(i+,-)), (Pi,Pi) G G 1 (p(i+, •)). 

We conclude 

JV(p(t+,-))=JV(p(t-,-)). 

The conclusion can be obtained in the same way if the wave is arriving to the junction 
from another road. Moreover, when two waves interact on the same road then there is a 
cancellation or gluing of waves and it is easy to check that N is constant or decreases. The 
proof is concluded. □ 
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6 Estimates on Flux Variation 

This Section is dedicated to the estimation of the total variation of the flux along a 
solution. We assume that every junction has at most two incoming roads and two outgoing 
ones. This hypothesis is crucial, because, as shown in the Appendix the presence of more 
complicate junctions provokes increase of the total variation of the flux. 

Lemma 6.1 Let f : [0, 1] — > M satisfy (J-l). Consider a network (Z,<J) in which each 
junction has at most two incoming road and two outgoing ones. Let p be a piecewise constant 
wave front tracking approximate solution. Then the map 

t>0^ Tot.Var.(f(p(t r ))), 

is not increasing. 



Proof. First of all we consider a single junction J with n < 2 roads with incoming traffic 
and m < 2 roads with outgoing traffic as in Section^. It suffices to study the case n = m = 2, 
the other ones are simpler. Let (pi.o, P4,o) be an equilibrium configuration in the junction 
J. Assume that a wave comes to the junction at the time t, we claim that 

Tot.Var.(/(p(i+, •))) = Tot.Var.(/(p(i- •))). (6.17) 

We begin assuming that the wave is on an incoming road, for example the first one, and that 
it is given by the values (pi, px t o). Let us define the incoming flux 

r(y)±\ f }«\ I*-*-** (6-18) 
I j(o-) if a < y < 1, 



and the outgoing flux 



I f{y) if o- < y < l. 



Clearly, since the wave on the first road has positive velocity, we have 

< Pi < a, f( Pl ) < f out ( Ph0 ). (6.20) 

Let (pi,...,p4,) the solution of the Riemann Problem in the junction J with initial data 
(pi, ...,P4,o) (see Theorem 3T). By definition (f(pifl),f(P2,o)) is the maximum of the map 
E on the domain 

^0 = {(71,72) € ^1,0 x ^2,o|^4 • (7i,72) T £ ^3,0 x ^4,0}, 
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and (/(/Si), /(/O2)) is the maximum of the map E on the domain 

H = I (71, 72) eiiix ^2 : o|^4 • (7i,72) T G ^3,0 x ^4,0}, 

where 

' [0,f"(^ )0 )] if i = 1,2, 



^7,0 — 



[0J out (p 3 fi)} if 3 = 3, 4, 



and ,by ( |O0| ), 

^i = [0,/ m (pi)] = [0,/( Pl )]. 

It is also clear that 

(/(Pl,0), /(P2, )) € ^0, (/(&,„), /(&,„)) G flft. 
To simplify the notations, define 

Oil = "31, "2 = "32 

then, by @, 

1 — ai = "41, 1 — a2 = "42 
We distinguish two cases. First we suppose that 

/(pi) < f(pifi), (6.21) 

(equality can not happen in the previous equation because the wave would have velocity 
zero). Then there results Cl C FIq, hence 

f(pi) < f(pi) < /(Pl.o), /(pi) + /(/5 2 )</(pi,o) + /(p 2 .o), (6.22) 

where the first inequality is due to the fact that the wave (pi , pi ) has negative velocity. We 
claim that 

/(p 2 ,o) < f{fr), (6.23) 

and 

/(/5 3 ) < fQ*,o), f(M < /(P4,o). (6.24) 

The points (f(pi,o),f(p2,o)), (f(f>i,o), f(f>2,o)) are on the boundaries of fio, fi respectively, 
where £7 is maximum, hence they are on one of the curves 

"171 + "272 = f° Ut (P3,o), (1 - "l)7l + (1 - "2)72 = f° Ut (p4fi), 72 = f m (p2,o)- 



Using ( p.21|) , we immediately get ( |6.23D . Let us assume that the two points are on the same 
curve, the general case being similar, for example on 

"171 + "272 = f° ut (p3,o)- 
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From (6.21) it follows that the map E is increasing on the curve 

7i ~ (71, /° ut (P3,o) - — 71), 

otherwise we contradict the maximality of E at (/(pi,o), /(p 2 ,o)) ■ Thus ai < a 2 , /5i = pi, 
and 

/(pi) = f(pi), /(AO = /(p 3 ,o) = / otlt (p 3 ,o). 

On the other hand, by fl6.22j ) and fl6.23| ), we have 

/(p 4 ) = (1 - ax)f{px) + (1 - a 2 )/(p 2 ) < 

< (1 - ai)(/(pi,o) + /(p 2 ,o) - /(&)) + (1 - a 2 )/(p 2 ) = 
= (1 - a 1 )(/(pi i0 ) + /(p 2) o)) + (ai - a 2 )/(p 2 ) < 

< (1 - ai)(/(pi,o) + /(p 2 ,o)) + («i - a 2 )/(p 2 ,o) = /(P4,o)- 
Using the Rankine Hogoniot Condition (2/7) in the junction ( |6.22| ), ( |6.23| ) and ( 6.24 ), we get 

Tot.Var.(/(p(t+,-))) = |/(pi)-/(p 1 )| + |/(p 2 )-/(p 2 ,o)| + |/(p 3 )-/(p3,o)| + |/(p4)-/(p4,o)| = 

= (/(p 2 ) - /(P2,0)) + (/(P3,0) - /(&)) + (/(P4,0) - /(At)) = 

= /(Pi.o) - /(pi) = /(pi,o) - /(Pi) = Tot.Var. (/(p(t- •)))• 
Suppose now that 

/(pi,o) < /(Pi), 
then f2o C and using the previous arguments 

/(Pl) = /(P1), /(P2) < /(p 2 , ), /(P3, ) < /(AO, /(P4 >0 ) < /(P4). 



By the Rankine Hogonot Condition in the junction (see (2.7)), we have 
Tot.Var.(/(p(t+,.))) = |/(p 1 )-/( / , 1 )| + |/(p 2 )-/(p 2) o)|+|/(p 3 )-/(p3,o)|+|/(p4)-/(p4,o)| = 

= (/(P2,0) - f(h)) + (/(AO ~ /(P3 )0 )) + (/(AO - /(P4,0))) = 

= (/(AO - /(pi,o)))+ - /(Pi) - /(pi,o) = Tot.Var. (/(p(i-, •))) . 

In the general case we have only to observe that the total variation of p does not increase 
on the roads (see || Chapter 6]) and when a wave approaches a junction we can use the 
previous argument, so the proof is concluded. □ 
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7 Solutions with a finite number of Big Waves 

In this section we prove existence and stability of solutions with a finite number of big 



Definition 7.1 We call T> n the set of all maps p : Uj/j i— > 1R defined on the network with 
bounded variation such that there exists a sequence {p v } of piecewise constant maps satisfying 



Tot.Var.(p v ) < Tot.Var.(p), 



N{p v ) < n, 



for each v € IN and 



Pv 



Notice that if p 6 T> n then 

Tot.Var.(p) < 



p vsxL . 



Tot.Var.(/(p)) +n. 



(7.25) 
(7.26) 

(7.27) 



The existence of solutions with values in the domain T> n is ensured by next Theorem. 

Theorem 7.1 Let f : [0, 1] — > 1R satisfying (^"1). Consider a road network (I,^J) in which 
all junctions have at most two incoming roads and two outgoing ones. Given n S IN and p € 



T> n , there exists an admissible solution p in the sense of Definition 2Jl such that p(0, •) = p, 
and p(t, ■) £ T> n , for each t > 0. 

Proof. Let {p^} C T> n be a sequence of piecewise constant maps such that 
Tot.Var.Q^) < Tot.Var.(p), p v — > p in L 1 , 



and p v the wave front tracking approximate solutions such that p u (0, •) = Pv- By Lemma 5.1 
we have 

MM e V n , t >0, v G JV 



and, by ( 7.27 ) and Lemma p.l[ 

Tot. Var.^ (*,•)) < 



Tot.Var.(/(^(t,-))) +n 



< 



< 



Tot.Var^/O^)) +n < 



||/'|U»Tot.Var.(p) 



Since the wave speeds are bounded, the maps p v are uniformly Lipschitz continuous from 
[0, +oo[ into L\ oc for the L 1 norm on every compact set, and are obviously uniformly bounded. 
Then by Helly Theorem (see || Theorem 2.4]), p v converge to some continuous map p E 
L ; 1 OC ([0, +oo[x Uj Ii,M) such that, up to redefining p on a set of zero measure, p(t, ■) has 
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bounded variation. Since for every t we can obtain p(i, •) as limits of p v {t ni -) for some 
t n — > t, we get p(t, •) 6 Z>„ for every f > 0. 

It is a standard argument, see ||, to prove that p solves the conservation law on the road 
network. Moreover, p(t, ■) G BV and one can easily check the other properties to guarantee 
that p is an admissible solution. The proof is concluded. □ 

Regarding stability of solutions we have: 

Theorem 7.2 Let f : [0, 1] — > M satisfy Consider a network in which all junctions 

have at most two incoming roads and two outgoing ones. Let n € IV, p and p be admissible 
solutions such that 

p(t,-),p{t,-) €V n nL\ t>0. 

There results 

\\ P (t,-)-p(t,-)\\ L1 < || P (o,-)-p(o,0IUi, 

for each t > 0. 

We begin proving a lemma. 

Lemma 7.1 Let f : [0, 1] — > 1R satisfy and J be a junction with two incoming roads 

and two outgoing ones as in Section |^. Let pi, pi.o, P4.o £ [0, 1], £ £ JR and x < fei fee swc/i 

- ^, c J [0,o-]\{pi,o} tf0<pi, <<7, 

x + 1 < oi, pi e < 

1 [0,T(p 1)0 )] */^<Pi,o<1. 



Lei p and p^ fee i/ie wave /ron< tracking approximate solutions of (1.1) on J such that 

pi(0, •) = pi • X]-oo,x] + Pi,o ■ X]x,bi]> Pi( Q > ■) = Pi ' X]-oo,2+e] + Pi,o ■ X]m+S,bi]> 
Pi (0,-)=pl(0,-)=p ii0 , i = 2,3,4. 

There results 

\\pf(t, ■) - Pl (t, < lei ■ Ipi.o - pi| = ||pf (o, •) - pi(o, -)]Ux 5 

/or eac/i t > 0. 

Proof. We begin assuming that (pi,o, ••• ) P4,o) is an equilibrium configuration. By possibly 
changing the notations, we can assume that £ > 0. Since we approximate the rarefaction 
fronts with many small shocks we have only to study the case in which the Riemann Problem 
(piiPifi) on the first road and the one (pi,o, ■■■,P4,o) m the junction generate only shocks. 
Let (pi, P4) be the solution of the Riemann Problem with initial data (pi, p2 ) o--., P4,o) (see 
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Theorem and Ai,o, Ai, A2, A3, A4 the velocities of the shocks generated by the Riemann 
Problems (pi,pi, ), (pi,Pi), (P2,p 2 ), {ps,p3,o), {pi, pi,o) respectively. 

If t > ^ ^ — - by the first part of the proof of Lemma 6.1 and the Rankine Hugoniot 

Ai,o 

Conditions on the roads (see Lemma 4.2]) we have 

\\ P {(t, •) - Pl (t, -)\\ L1 = J*L (\ P1 -h\- |Ax| + E k° - ^1 • M 



,^(l/(pi) ~ f(Pi)\ + E l/(^.o) - f(Pi)\) = ATot.Var.(/(p(t, •))) = 
k1 r Tot.Var.(/(p(0,-))) = A|/G>i)-/G>i,o)| = |£|-|pi-pi,o| = \\A (0, -)-Pi(0, OIU- 



|Ai,o| ' ' |Ai, o| 

In the case in which (pi^o, ■■•,P4,o) is not an equilibrium configuration we have only to 
recall that the L 1 — distance between the solutions decreases on each roads (see of |], Corollary 
6.1]) and use the same arguments. 

This concludes the proof. □ 



Proof of Theorem 7.2. Let p and p be front tracking approximate solutions such that 

11/(0, -)-p(0, Oik <~, \\p k (0,-)-p(0,-)\W<~, keN, 

Tot.Var.(p fe (0, •)) < Tot.Var.(p(0, •)), Tot.Var.Q5 fe (0, •)) < Tot.Var.(p(0, •)), k e IN. 
Now consider finitely many wave front tracking approximate solutions p k '°, ...,p ,N 

k,0 _ k k,N _ ~k 

P = p , p = p 

where p fc '' 1 is obtained by p k ' h ~ 1 shifting and rescaling only one jump as in M and ||. 
Precisely denoting 

AT-l 

1=0 

there exist \,£ £ 1R and fe {1, ...,iV — 1} such that 

1-2 N+l 

with 

ft-ft-i<Z<ft+x-fa 
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for each i = 1, N. In this way we have 

N 

£ 11^(0, •) - p fe >' l+1 (0, -)\\ L1 = \\ P k (0, .) - p*( , -)IUx- (7.28) 
fc=i 

Since the distance between solutions decreases on each road (see ||, Corollary 6.1]) and by 
the previous lemma, we have 

\\ P k > h (t, •) - P k > h+1 (t, .)|| £ i < ||p^(o, •) - p^ i+1 (o, oiUi. 

So, by 

■) - ~ P \t, oiUi < ^ •) - p fc - h+1 (t ) 011^ < 

iV 

< E 11^' •) - 1+1 (°> OIU- = ll/(o, •) - P fc (o, OiUi. 

Moreover there exists a decreasing sequence {k n } C iV such that p' c " — > p and /5 fc ™ — > p in 
L 1 as fc n — > +oo. Hence 

|| P (t, -) - -)IUx < ]| P (o, -) - p(o, -)IU-, 

as to be proved. □ 

8 Existence and Stability of Solutions in L 1 

Let us first consider the case in which (Tl) holds true. 

Theorem 8.1 Let f ; [0, 1] — > 1R satisfy (Fl). Consider a road network in which each 
junction has at most two incoming roads and two outgoing ones. Let p be an initial data 
in L\ oc and fix T > 0. Then there exists a unique admissible solution p defined on [0,T], 
obtained as limit of wave front tracking approximate solutions such that p(0, ■) = p. Moreover 
if p G L 1 then p{t, ■) G L 1 . 

If p and p are admissible solutions obtained as limit of wave front tracking approximate 
solutions such that p(t, ■), p(t, ■) G L , for every t > 0, then for each t > 

||p(V)-p(V)IUi < ||p(0, -) - p(0, (8.29) 

Proof. We begin proving the existence of a solution for p G L 1 . There exists {p n } sequence 
of piecewise constant maps defined on the network such that 

p n eV n HL\ pn—>P in L 1 . (8.30) 
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Let p n be the wave front tracking approximate solutions with p n (0, •) = p n . Fix n < m, by 



Lemma 5.1, there results ■) € 2?n C T> m ,p m (t, ■) G 2? m and by Theorem 7.2 



Hence {p n (t, •)} is a Cauchy sequence in L ([0,T] x Ft, Ft). Then there exists p such that 
Pn(t, — * hi L 1 ([0,T] x Ft, Ft). It is easy to check that p is as admissible solution. 
The case of L\ oc can be obtained by localization. 

Now we prove ( |B.29 ). Let {p«}, {p n } be sequences of wave front tracking approximate 
solutions such that p n (t, •), p n (t, ■) G P n H L 1 and p„ — ► p, p n — > p in L 1 ([0,T] x Ft, Ft). 
By Theorem 7.2, we have 

IM*, ') " Pn(t, < l|Pn(0, •) " p„(0, -)|Ul. 



Therefore J8.29|) is proved and uniqueness holds true. 

We now relax the assumption (Tl), namely we suppose that / satisfies (F). 
Let {fv} be a sequence of maps satisfying (J-l) such that 



and 



fu 



f v {p) = max Ju(p), v G IV 
P e[o,i] 



/ in L°°([0, 1]) and jj — /' in L°°([0, 1]). 



□ 



.31) 



.32) 



Moreover let p be an initial data in Lj oc . We know that there exists a unique p v = p u (t,x) 
admissible solution to the Cauchy Problem on the network (see Theorem 8.1) obtained as 
limit of front tracking approximate solutions for 



Pt + fu(p)x = 0, 



v G IN. 



.33) 



Theorem 8.2 Let f : [0, 1] — > M satisfy (J-). Consider a road network in which all the 
junction have at most two incoming roads and two outgoing ones. Let p be an initial data in 
L\ oc and fix T > 0. Then there exists a unique admissible solution p defined on [0, T], with 



p(0, •) = p, obtained as limit of solutions to ( 8.35 ). The limit does not depend on the choice 



of the functions /„ and if p G L 1 then p(t, •) G L 1 . 

If p and p are such admissible solutions and satisfy p(t, ■), p(t, •) G L , for every t > 0, 
then 

||p(*,0-p(V)IU* < ||p(0,-)-p(0,-)|| L1 . (8.34) 



The Theorem can be proved exactly as Theorem 8.1 from next Lemmas 
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Lemma 8.1 Let f : [0, 1] — * M satisfy (J-). Consider a road network in which all junctions 
have at most two incoming roads and two outgoing ones. Let p be an initial data in T> n n L 1 
and fix T > 0. Then there exists a unique admissible solution p defined on [0,T], with 
p(0, •) = p, obtained as limit of solutions to ( 8.3$ ). The limit does not depend on the choice 
of the functions /„ and and p(t, ■) £ L 1 . 



Proof. Let {/„} be a sequence of maps satisfying (J-l), ( |8.31| ) and (3.32) and p v be the 



admissible solutions for the Cauchy problems associated to /„. By Theorem 7.1 we have 

p v (t, •) e V n , veN, 0<t<T. (8.35) 

Moreover there results 

\\p v (t, •) - p^t, -)\\ L t <C-\\fJ - / M '|| L « • Tot.Var.(p), (8.36) 



where C depends only on /. By ( fe.32j ), {p v } is a Cauchy sequence in L 1 , then there exists 
p G L 1 such that p v — > p in L 1 . Moreover, p is an admissible solution and satisfies p(0, •) = p. 



From (8.36) we have that p does not depend on the choice of {f„}, so we are done. 



□ 



Lemma 8.2 Let f : [0, 1] — ► M satisfy Consider a road network in which all junctions 

have at most two incoming roads and two outgoing ones. Fix T > and let p, p be admissible 
solutions in L 1 , obtained as limit of solutions to /[8.3$) , defined on [0, T]. If p(0,-), p(0, ■) G 
T> n , then 

\\p(t,-)-p(t,-)\\ Ll < ||p(0,.)-p(0,.)|| £ i, 

for each < t < T . 



Proof. Let {f v } be a sequence of maps satisfying (8.31) and (8.32) and p v be the 

admissible solutions associated to /„ such that 



p v (0,-)=p(0,-), p u (0,-) = p(0,-). 



.37) 



By Lemma iA, we have p v — > p and p v — > p in L 1 . By Theorem 8^ and ( 8.37 ), for each 
< t < T and v G IV, there results 

IIp(*, •) - OIU 1 < IIp(*, - ^(*: OIU 1 + •) - OIU 1 + l!- 5 ^ ') ~ Mt, -)\\ L i < 

< \\p(t, ■) - p v {t, .)|| L i + 11/9(0, •) - p(0, OHli + •) - P,(t, .)\\ L1 — . 11/9(0, •) - p(0, OHli. 
So the proof is concluded. □ 
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9 Time Dependent Traffic 

In this section we consider a model of traffic including cross lights and time dependent 
traffic. The latter means that the choice of drivers at junctions depends on the period of the 
day, so during the morning the traffic flows towards some direction and during afternoon it 
may change towards another direction. This means that the matrix A depends on time t (see 
Section |). 

Consider a single junction as in Section || with two incoming roads and two outgoing ones. 
Let ot\ = a±(t), a 2 = a 2 (t) be two piecewise constant periodic functions such that 

a 1 (t)^a 2 {t), (9.38) 

for each t > 0. Moreover let xi — Xi Wj X2 — Xzif) be piecewise constant periodic maps such 
that 

Xi(*)+Xa(*) = l, Xi(*)e {0,1}, i = l,2 

for each t > 0. The two maps represent traffic lights, the value corresponding to red light 
and the value 1 to green light. 

Definition 9.1 Consider p = p(t, xi, X4) = (pi(t, x±), p^{t, £4)) with bounded varia- 
tion. We say that p is a solution at the junction J if it satisfies (i), (ii), (iv) of Definition 
\2.\ and the following property holds: 

(v) f(p 3 (t,a 3 +))=a 1 (t)xi(t)f(pi(t,b 1 -))+a 2 {t)x2{t)f(p2(t,b 2 -)) and 

f(p 4 (t,a 4 +)) = (1 - a 1 (t)) X i(t)f(pi(t,b 1 +)) + (1 - a 2 (t)) X 2(t)f(p2(t,b 2 +)) for each 
t > 0. 

Assume that at time t one of the maps cxi(-), a 2 (-), Xi('); X2(') jumps, then we have to 
solve a new Riemann Problem in the junction hence four waves are generated and 

N(p(i+,-)) <N(p(t-,-)), (9.39) 

Tot.Var.(/(p(f+, •))) < Tot.Var.(/(p(t- •))) + 4/(a). (9.40) 
Then the map N(p(t,-)) is still non increasing while 

Tot.Var.(/(p(i 2 ,-))) < Tot.Var.^^, •))) +4f(v)$(t 1 ,t 2 ), 

for each < t\ < t 2 , where 

®(h,t 2 ) = ^ <tl< <2, X% jumps in ti} + #{ti\ti < t t < t 2 , oti jumps in ti}^. 

i=l,2 
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Therefore, for fixed T > 0, we have uniform bounds of the total variation on the interval 
[0, T], and using arguments as in the previous sections we obtain existence and stability for 
the Cauchy Problem. However, the total variation of f(p) does not depend continuously 
on the total variation of the maps cxi(-), o^G)- Indeed consider a single junction with two 
incoming roads and two outgoing ones without traffic lights, i.e. \i = 1, an d let 



where < (5% < Pi < \ and < t < T. Consider the initial data (pi,o, P2,o, P3,o, Pi.o), where 



This is an equilibrium configuration for the choice on = Pi, i = 1, 2, hence the solution of the 
Riemann Problem is identically equal to the initial data for < t < t. At time t = i we have 
to solve a new Riemann Problem. Let (pi, p2, p3, Pa) its solution, there results (see figure 2) 




/(Pl,o) = /(P4,o) = /(O"), 



/(P2,o) = /(P3,o) 



1-P2 



fifh) = f(P4) = f{&) 



f(pl) = f(h) = 



I -Pi 



f(o ) 




(i-p I )Y + d-P 2 ) y= f(o ) 




f(o ) 



Y 



figure 2 



Now, if Pi — > /3 2 , then 



Tot.Var.(ai; [0,T]) — > 0, Tot.Var.(a 2 ; [0,T]) — > 0, 



but 




hence Tot.Var.(/(p); [0, T]) is bounded away from zero. 
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A Appendix: Total Variation of the Fluxes 

In this section we show an example in which the total variation of the flux increases due 
to interactions of waves with junctions. 

Consider a single junction with three incoming roads and three outgoing ones, the matrix 



A = 



/I I I \ 

/ 2 2 3 \ 
111 

3 2 2 

u o y 




(A.41) 



figure 3 

and constants pi, pi.o, ■ p6,o € [0, 1] such that 

Pi,o = P3,o = P4,0 = P5,0 = cr, a < p 2 ,o < 1, < p 6) o, Pi < cr, /(p 2 ,o) = ^, /(P6,o) = g- 

Assume that /(a) = 1, then (pi.o, P6,o) is an equilibrium configuration and p given by 

, n v / Pi,o if a;i < a; < 6i, . 
Pi(0,a;) = ^ pi(0, •) =Pi,o, « = 2,..., 6, 

I pi if x < Xi, 



is a solution. Moreover the plane 



6 71 + 6 73 = 1 



does not intersect the cube [0, l] 3 and the point (f(pifl), /(pe,o)) is on the intersection of 
the planes 

111 111 

27i + 2T2+g73 = l, 371+ 272+ 273 = 1, 



that is the line described by the map 

7i ^ (71,2 - 3-71,71)- 



(A.42) 



At some time say t the wave (pi,pi,o) interacts with the junction. Let (pi, ....,pe) be the 
solution of the Riemann Problem at the junction for the data (pi, p2,o, P6,o)- Since the 
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map E increases on the line described by (A. 42), the point (f(pi), f{p&j) is on the curve 



(A.42) and 



f(jk) = m) = fo>i), fih) = 2 - §/(pi), am = m) = f(°), f(M = lf(pi)- 



We get 



while 



Tot.Var.(/(p(i-,-))) =l-/(pi), 
Tot.Var. (/(p(t+, •))) = 4(l - f( Pl )) > Tot.Var. (/(p(i- •)))• 



B Appendix: Total Variation of the Densities 

Consider a junction J with two incoming roads and two outgoing ones that we param- 
eterize with the intervals ] — oo, b\], ] — oo, 62], [03, +oo[, [04, +oo[ respectively. Fix the 
constants a±, ot2 such that < ol\ < 0L2 < 1/2 and set 

0:3,1 = OLi, 013,2—0:2, Q!4i = l — ai, 013,2 = l-a 2 . 




Define a solution p by 

pi if xi < a; < 61, 



pi{0,x) 



pi ii x < %i, 



figure 4 



p 2 (0,x) = p 2:0 , P3(0,x) = p 3:0 , P4(0,x) = p 4 , , 

(B.43) 



where pi, pi,o, p2,o, P3,o, P4,o are constants such that 

O < P2,0- P3,0 < 1, < pi < a, pi, = P4,0 = c, 

01 



(B.44) 



/(Pl,o) = /(/>4,o) = /(c), /(P2,o) = /(P3,o) 



1 — «2 



/(c), 
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so (pifl, P2,q, P3,0i P4,o) is an equilibrium configuration. 

After some time the wave (pi,/?i : o) interacts with the junction. Let [pi, p2, P3, pi) be 
the solution of the Ricmann Problem in the junction for the data [p±, P2,0j P3,0i P4,o)- By 
flB.43|) and (|B.44D there results 



/(pi) = /(pi), /(AO 



1 — «2 



/(Pa) ^ r^/C* 



cv 2 



and 



1 — a 2 1 — «2 

< p 3 < ct < /5 2 < 1 



(B.45) 




figure 5 



Therefore if pi — > pi_o = cr then 

/(Pa) 



(Vl 



1 — «2 



/(a) = /(p 3 , ), 



and by (B.45), (B.44) we have p 3 — > t(p 3 ). Therefore we are able to create on the third 
road a wave with strength bounded away from zero using an arbitrarily small wave on the 
first one. 
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